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Abstract 


A nonperturbative analytic solution of the high charge and energy (HZE) 
Green’s function is used to implement a computer code for laboratory ion beam 
transport. The code is established to operate on the Langley Research Center 
nuclear fragmentation model used in engineering applications. Computational 
procedures are established to generate linear energy transfer (LET) distributions 
for a specified ion beam and target for comparison with experimental measure- 
ments. The code is highly efficient and compares well with the perturbation 
approximations . 


Introduction 


by a mean energy E 0 and energy spread cr such that 


Green’s function has been identified as the likely 
means of generating efficient high charge and energy 
(HZE) shielding codes for space engineering which 
are capable of being validated in laboratory experi- 
ments (ref. 1). A derivation of Green’s function as a 
perturbation series was promising for development of 
a laboratory- validated engineering code (ref. 2), but 
computational inefficiency provided a major obstacle 
to code development (ref. 3). More recently, non- 
perturbative approximations to the HZE Green’s 
function have shown promise in providing an effi- 
cient validated engineering code (ref. 4). Described 
herein is a laboratory code that uses a nonperturba- 
tion Green’s function to derive LET spectra for ion 
beams with an atomic number Z < 28. These ions 
correspond to the major components of the galactic 
cosmic ray spectrum. 


W E) = VZb cxp 

The usual method of solution is to proceed in solv- 
ing equation (1) as a perturbation series (refs. 1 
and 5). In practice the computational requirements 
limit the usefulness of the technique for deep pene- 
tration (ref. 3). 

The Green’s function is introduced as a solution 
of 
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Theory 


subject to the boundary condition 


We restrict our attention to the multiple charged 
ions for which the Boltzmann equation may be re- 
duced to 
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( 1 ) 

where (f)j(x>E) is the flux of ion type j at x with 
energy E (in MeV/amu), Sj(E) is the change in E 
per unit distance, a j is the total macroscopic reaction 
cross section, and aj^ is the macroscopic cross section 
for collision of ion type k to produce an ion of type j. 
(See refs. 1 to 5.) The solution to equation (1) is 
found subject to the boundary condition 


,E) = fj(E) (2) 

which for laboratory beams has only one value of j 
for which fj(E) is not zero, and fj{E) is described 


Gj m {0, E : E 0 ) ~ 6 jm 6(E - E 0 ) (5) 

where Sj m is Kronecker’s delta and <5(J57 — E 0 ) is 
Dirac’s delta. The solution to equation (1) is given 
by superposition as 

4>j(x, E) = Y2 [ G Jk(x, E, E’)f k {E') dE' (6) 

k J 

If Gjk(X) E, E f ) is knowm as an algebraic quantity, 
then equation (6) may be evaluated by simple inte- 
gration techniques and the associated errors in nu- 
merically solving equation (1) are avoided (ref. 6). 

The above equations can be simplified by trans- 
forming the energy into the residual range as 

J Jo Sj(E') 


( 7 ) 


and defining new field functions as 


In the above equation, 


ipj(x,rj) = Sj(E)<j>j(x,E) 

Q jm (x, tj , r m ) = Sj[E)Gj m (x, E, E ) 
fj( r j) = Sj(E)fj(E) 


Equation (4) becomes 
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dx drj 
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with the boundary condition 

&jm{Q‘> r j'> r m) = ~ r m ) 

and with the solution to the ion fields given by 


( 12 ) 


4>j{x,rj) = E / Sjrn(x,r J ,r , m )f m {r , m ) dr' m (13) 
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Note that vj is the range scale factor such that vjtj = 
i/ m r m and can be expressed as vj = Zj/Aj. The 
solution to equation (11) is written as a perturbation 
series: 


g(j) = exp {-(Jjx) (19) 


and 
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Note that <7^(x, 5 r' R ) is purely dependent on x for 

u\ 

i > 0, and we represent this expression as QF{x). 
(See ref. 3.) In terras of equations (14) to (20), the 
solution to equation (1) becomes (ref. 3) 

ipj(x,rj) = exp (~<Tjx) fj(rj + x) 

+ - Sn( r m,tt)] 

m } i 

(21) 

In equation (21), r[ n u and r* t are given by the 
upper and lower limits of equation (17). The symbol 
F m (r f m ) refers to the integral spectrum 

Tm( r m) = [. fm( r ) dr (22) 

JrL 


Qjrnix, Tj,r' m ) = &jm( X ’ r j’ r m) ( 14 ) 
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where 

gf n {x,r 3 y m ) = g{j)6 jm S(x + Tj - r' n ) (15) 


and 
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where 



is zero unless 


(16) 


— (rj + x) < r' m < Tj + x (17) 

V m Vjji 

for Vjji > vj. If Vj > v m , as can happen in neutron 
removal, the negative of equation (16) is used and the 
upper and lower limits of equation (17) are switched. 
The higher terms are approximated as 
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- E 
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We note that 


F m (r m ) = Fm(E') 


with 


and 


roc 

F m (E') = / fm{E) dE 
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We now introduce nonperturbation terms for the 
summation in equation (21). First, we recall that the 
g - function of n arguments is generated by the pertur- 
bation solution of the transport equation neglecting 
ionization energy loss (ref. 1) given by 

^ ^ s &ik 9km{ (26) 

subject to the boundary condition 

£jm(0) — (27) 
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The solution is 

9jm{ x ) ~ ^jm#( m ) T O'jmdU't rn ) 

+ X a 3k <r km9(j , k,m) + ■■■ (28) 

k 

It is also true that 

9jm(x) = '529jk(x - y)9km(y) ( 29 ) 

k 

for any positive values of x and y. Equation (29) 
may be used to propagate the function gj m (x ) over 
the solution space, after which 

Gjm(x, r ji r m) ~ exp(-o-jx) 6 jm 6(x + rj - r' m ) 

, Vj [< q jm (x ) - cxp(-(7jx)(5j m ] 

• / \ 

X\Vm ~ Vj) 

The approximate solution of equation (1) is then 
given by 


ipj(x,rj) = exp (-<Tjx) f (rj + x) 


+ E 


femW - exp(-(T ? x) 6 jm 
%(Vm. k'j) 



which is a relatively simple quantity (ref. 4). 


(31) 


Comparison With Perturbation Theory 

The first step in testing the new code is compar- 
ing its results with the perturbation theory results. 
The perturbation code has been previously compared 
with 20 Ne transport experiments at the Lawrence 
Berkeley Laboratory (LBL) BEVALAC accelerator 
(refs. 7 and 8). The comparison required the use 
of acceptance functions to obtain detector response, 
and results of those comparisons are shown in fig- 
ure 1. The dynamic range of the detector system was 
inadequate for the Be and B ions especially at the 
lower depths. The results are thought to be within 
30 percent over the dynamic range of ions for which 
the detectors were designed (ref. 8). The primary er- 
rors in the computation are attributed to the nuclear 
cross sections and approximations used in applying 
the acceptance functions. A sample LET spectrum 
for the C ions is shown in figure 2, which indicates 
excellent agreement between theory and experiment. 
The spectrum for the Ne ions is shown in figure 3. 
The agreement is less favorable since Ne fragment 


isotopes have not been evaluated in the computa- 
tions. (See ref. 8 for Ne fragment isotope effects 
at 15 cm.) The main uncertainty is believed to be 
the fragmentation cross sections. The only complete 
sets of fragmentation cross sections measured are for 
12 C and 1? 0 projectiles at both 1.05 GeV/amu and 
2.1 GeV/amu, and 20 Ne cross sections are obtained 
(ref. 9) as an extrapolation. The Ne, F, and N cross 
sections scaled according to reference 9 are among 
the more reliable unmeasured cross sections. 

We converted the perturbation code to use the 
NUCFRAG data base (refs. 10 and 11). Direct 
comparison of the perturbation code and the non- 
perturbation Green’s function code is then possible. 
Figure 4 shows the first collision term from both the- 
ories for a 600 MeV/amu 56 Fe beam at three depths 
in water. The differences in spectral shape are due 
to simplification of the attenuation term in the non- 
perturbation theory (ref. 1). The total ion flux of 
each type is the same as we have shown elsewhere 
(refs. 1 and 6). The nonperturbation term represents 
the average spectrum, while perturbation theory re- 
tains the spectral shape. The second collision terms 
are shown in figure 5. There is better spectral de- 
tail maintained in the perturbation code than in the 
nonperturbation code, but the same total ion flux 
for each ion type is obtained for both. The (approxi- 
mate) third collision term of the perturbation theory 
is shown with the remaining nonperturbation terms 
in figure 6. The total ion fluxes for the two theo- 
ries are nearly equal for heavy fragments but differ 
substantially for lighter fragments. By comparing 
figures 4, 5, and 6 we see that the sequence of per- 
turbation terms appears to be converging to a re- 
sult similar to the nonperturbation result. Although 
perturbation theory converges rapidly for the heavier 
fragments, the light fragment production from higher 
order terms is important for this heavy ion beam and 
cannot be neglected. This deficiency in accounting 
for light fragments is not a great concern since these 
terms are efficiently evaluated by the present non- 
perturbative techniques. Results for the first three 
terms of perturbation theory are compared with re- 
sults for all terms of the nonperturbation theory in 
figure 7. The nonconvergcnce of the light fragments 
in the first three perturbation terms is most apparent 
by comparing results for the third-order and higher 
nonperturbation terms (fig. 6(c)) with results for the 
third-order perturbation term (fig. 6(d)), for exam- 
ple. Aside from the convergence issue, the compu- 
tation time required for the nonperturbation code is 
approximately 10 minutes on a VAX 4000/500 sys- 
tem, compared with 15 minutes for evaluation of 
the first collision term, 45 minutes for the second 
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collision term, and 90 minutes for the third colli- 
sion term for the perturbation code. We felt that 
little would be gained by evaluating the fourth term, 
which would require more computer time than the 
first three terms combined. In addition to the ^Fe 
ion beam results given in figures 4 to 7, the nonpcr- 
turbative method generates the Green’s function for 
any ion at Z < 28 for any arbitrary initial energy. 
Each additional ion or energy requires an additional 
150 minutes to evaluate the first three collision terms 
using the perturbation-based code. This should am- 
ply demonstrate the power of the nonperturbative 
solutions. 


A Practical Green’s Function Code 

Although the nonperturbation Green’s function 
provides a rapid computational method, the spectral 
terms are replaced by averages over the spectral do- 
mains of higher order terms. These average spectra 
can be easily corrected with the spectral distributions 
using perturbation theory. The first collision Green’s 
function is given (ref. 2) as 


Gj 1 J l ( x ’ r J’ r m) = Z r ,,"i exp(-<7jX, - a m x m ) 

1 Vm ~ ' Vj | 

(32) 

We rewrite Q ^ in terms of its average value (given 
by eq. (16)) as 


+ b jni( x )( T m ~ f 'm ) ( 33 ) 


U in - Vi X 


where 


_/ r m/ r m,u 


is the midpoint f' m between its limits given by equa- 
tion (17). The bj m term of equation (33) has the 
property that 



bjm{x)(r' - r' m ) dr' = 0 


(34) 


to ensure that the first term of equation (33) is indeed 
the average spectrum as required. The spectral slope 
parameter is 


bjm{ x ) — 


m [exp ( &m x ) Cxp( 


x{v m 


Vj)\v m - Vj [ 


(35) 


the higher order terms. The spectral correction given 
in equation (33) is included in the present Green’s 
function code. 


LET Spectra for Laboratory Beams 

We use the boundary condition appropriate for 
laboratory beams given by equation (3). The cumu- 
lative spectrum is given by 


Fj(E) 



( 36 ) 


The cumulative energy moment needed to evaluate 
the spectral correction is 


Ej(E) = -E 0 


1 — erf 


( E-E c 

V V2a 


+ 


v/27r 


exp 


[E-E 0 

2a 2 


(37) 


The average energy on any subinterval (E \ , E 2 ) is 
then _ __ 


EjjEj) - EjjEt) 

Fj(E\) - Fj(E 2 ) 


(38) 


The beam-generated flux is 


= cxp(-ajx) }j{r j + x) 

m,i 

+ X/2 (l) Vm(E) -4] pm(4,u) - £»(4,,)] 

m 

(39) 

where E is evaluated using equation (38) with E\ 
and E 2 as the lower and upper limits associated with 

r 'mJ and r 'm,u- 

The differential fluencc spectra for a 600 MeV/ 
amu 56 Fe beam with a 2.5 MeV/amu standard devi- 
ation are shown in figure 9. The results in figure 9 
are similar to those shown in figure 8 except the ends 
of the spectra arc rounded by the energy spread of 
the primary beam. The LET distribution is found 
using the methods of reference 12. The correspond- 
ing LET spectra are shown in figure 10. The highest 
LET peak is due to the primary beam and the iron 
fragments. The successive peaks below iron are due 
to lower atomic weight fragments. 


Results for the nonperturbation Green’s function 
with spectral corrections are shown in figure 8 and 
should be compared with the results in figure 7. A 
similarly simple spectral correction could be made to 


Concluding Remarks 

The present analysis has resulted in a new 
transport code for high charge and energy (HZE) 
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ions which is a highly efficient nonperturbative an- 
alytic solution. Indeed, the basic solution for 
every ion through 58 Ni is generated in minutes 
on a VAX 4000 series computer. The solutions 
compare favorably with the perturbation series so- 
lutions, which have been validated to some ex- 
tent in laboratory experiments with 20 Ne beams. 
The code holds promise of an efficient space en- 
gineering code that can be tested in laboratory 
experiments. 

NASA Langley Research Center 
Hampton, VA 23681-0001 
June 23, 1993 
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— Numerical code (ref. 8) 

Perturbation theory (1st generation) 
Perturbation theory (2nd generation) 
* Experiment 




Figure 1. Integral fluence for all fragments produced by neon nuclei incident on various thicknesses of water 
as function of water column thickness. (From ref. 6.) 


6 


— Perturbation theory (2nd generation) 
Perturbation theory (1st generation) 
• Experiment 



Figure 2. Differential fluence for carbon nuclei produced by neon nuclei incident on various thicknesses of water. 
(From ref. 6.) 
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Figure 3. Differential fluence for neon nuclei incident on various thicknesses of water. (From ref. 6.) 
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(a) Nonperturbation theory; water depth, 5 cm. 



(b) Perturbation theory; water depth, 5 cm. 


Figure 4. Differential fluence (first collision) for 600 MeV/amu 56 Fe beam incident on water slab from 
nonperturbation and perturbation theories. 
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(f) Perturbation theory; water depth, 15 cm. 
Figure 4. Concluded. 
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Figure 5. Differential fluence (second collision) for 600 MeV/amu 56 Fe beam incident on water slab from 
nonperturbation and perturbation theories. 







(a) Nonperturbation theory; water depth, 5 cm. 



(b) Perturbation theory; water depth, 5 cm. 


Figure 6. Differential fluence (higher order collision) for 600 MeV/amu 56 Fe beam incident on water slab from 
nonperturbation (all terms) and perturbation (third term only) theories. 
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(d) Perturbation theory; water depth, 10 cm. 
Figure 6. Continued. 





(a) Nonperturbation theory; water depth, 5 cm. 



(b) Perturbation theory; water depth, 5 cm. 


Figure 7. Differential fluence for 600 MeV/amu 56 Fe beam incident on water slab from nonperturbation (all 
terms) and perturbation (first three terms) theories. 
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(c) Nonperturbation theory; water depth, 15 cm. 



(f) Perturbation theory; water depth, 15 cm. 
Figure 7. Concluded. 
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(a) Water depth, 5 cm. 



(b) Water depth, 10 cm. 


Figure 8. Differential fluence for 600 MeV/amu 56 Fe beam incident on water slab from spectral corrected 
nonperturbation theory. 





(c) Water depth, 15 cm. 
Figure 9. Concluded. 
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(b) Water depth, 10 cm. 

Figure 10. LET distribution for differential fluence for 600 MeV/amu o6 Fe beam with 2.5 MeV/amu standard 
deviation incident on water slab. 
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